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Charge neutral bilayer graphene has a gapped ground state as transport experiments demonstrate. 
One of the plausible such ground states is layered antiferromagnetic spin density wave (LAP) state, 
where the spins in top and bottom layers have same magnitude with opposite directions. We propose 
that lightly charged bilayer graphene in an electric field perpendicular to the graphene plane may 
be a half metal as a consequence of the inversion and particle-hole symmetry broken in the LAF 
state. We show this explicitly by using a mean field theory on a 2-layer Hubbard model for the 
bilayer graphene. 
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Half metals are a class of materials in which electrons 
with one spin orientation are metallic and electrons with 
opposite spin orientation are insulating [HQ- In a half 
metal, electrical current can be fully spin polarized. This 
property is attractive in spintronics [3!-|5|. Possibility of 
graphene-based half metal has been interesting for its po- 
tential application in electronic devices. Soon after the 
discovery of graphene. Son, Cohen and Louie have ap- 
plied first principles calculations to propose half metal- 
lic phase in a zigzag graphene nanoribbon with external 
transverse electric fields |6J] . Due to the technical difficul- 
ties in the devices, the predicted half metal in graphene 
ribbon has not been confirmed yet in experiments. Here 
we propose that bilayer graphene (BLG) in an electric 
field perpendicular to the graphene plane can be a half 
metal. Our proposal is based on the recent works to sug- 
gest that the ground state of BLG at half filled (charge 
neutral case) may be a layered antiferromagnetic spin 
density wave (LAF) state with a gap of about 2 meV. 
In the half filled LAF state, the spins are polarized op- 
positely on the top and bottom layers, which break the 
spin SU(2) symmetry, but are invariant under a com- 
bined transformation of time reversal (T) and inversion 
(/), as we can see in Fig. 1(a). The LAF state also con- 
serves the combined symmetry of inversion and particle- 
hole conjugation (C), where the operation of C changes 
an electron operator to a hole on the same site, which 
also reverses the sign of the magnetic moment of a par- 
ticle (or hole). The electric field breaks / and C. Their 
combination I ® C is conserved in the neutral graphene, 
but is broken in the charged graphene, under the electric 
field. Therefore, an electric field applied to the charged 
graphene opens the possibility for a half metallic phase 
with a net magnetization. We examine the half metallic 
BLG by applying a mean field theory on a 2-layered Hub- 
bard model. Our prediction may be tested in the BLG 
device with double gates J]]. 

We start with a brief summary of the recent works 



on the BLG|7h37I|. Theoretically, in the single electron 
band picture the BLG is a gapless semiconductor with 
parabolic valence and conduction bands touching at the 
high symmetry points K and K' . The state is unstable 
in the presence of electron interaction. Experimentally, 
there are clear evidences that the BLG at half filled has 
a gapped ground state f?!. [l6l - Eoj |. Velasco et al. |7| have 
applied a perpendicular electric field on a high quality 
suspended BLG. The energy gap is found to decrease as 
the field increases and to close at a field lbmV'nm~^ . 
Two of the most promising states consistent with the 
gapped ground states measured in transport experiments 
are the LAF [H, ll?!, 'ss', 'iol, 'il'l and the quantum spin 



Hall states 26|, |27|, |33, 38, 39]. It will be important and 
interesting to explore the possible experimental conse- 
quences of these states and determine the true ground 
state of the BLG. In the present Letter, we predict a half 
metallic phase in the LAF state of the BLG in the pres- 
ence of an electric field at slightly charged graphene. Our 
result may be used to resolve the controversial issue of 
the ground state in BLG, and should be of importance to 
the potential application of the graphene in spintronics. 

Wc consider a BLG system in an applied perpendicular 
electric field. The Hamiltonian is given by 



H — Ho + Hi, 



H„ 



(1) 



where Hq = H intra + Hinter IS thc kinetic energy part, 
Hu is the on-site Coulomb interaction, and Hp describes 
the electric potential due to the applied electric field. The 
intra-layer hopping is given by 



Hi, 



l{ij)a- 



fJ-^ni„{i) (2) 



Here, aia (bia-) are the electron annihilation operator on 
sublattice A (B), and ^ = 1,2 for the bottom and top 
layers as illustrated in Fig. 1, a =t,i and i (j) denote 
the spin and site, respectively. () sums over all nearest 



neighbor sites. Here we only consider the nearest neigh- 
bor hopping for simphcity and expect the small remote 
hopping will not change the basic physics, fx is the chem- 
ical potential, /i = corresponds to the half filling, or the 
charge neutral point. At /i > 0, the chemical potential 
acroses the conduction band. We consider an interlayer 
hopping between two sites i and i' on-top to each other, 
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[b\^{i)a2a{i') + h.c. 



(3) 



The Hubbard U term is given by Hu — 
UY^i^ni-^[i)ni]^{i). The effect of the external elec- 
tric field Eq is modeled by an electric potential V 
between the two layers, 



Hp = y^yini„{i), 



(4) 
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with Vi = {—iyV/2, and V is related to Eq as below. 
Note that the electron charge on the two layers may 
be redistributed in the presence of i^oO, lil- Let E 
be the electric field between the two layers, and assume 
the graphene sheets to be infinitely large 2-dimensional 
planes, then we have V — +edoE, with — e the electron 
charge, do the interlayer distance. E is related to Eq by 



E = Eq - 27re(p2 - Pi), 



(5) 



where the electron density in layer / is given by pi = 
(Scr i nicrii)) / S , with S the area of each layer, and (Q) is 
the average value of operator Q. In the above equation, 
we have assumed the dielectric constant for the BLG to 
be 1 as suggested in the literaturelillii]. 

We use a mean field approximation for the Hubbard 
term, and solve the Hamiltonian self-consistently. 



(6) 
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There are four atoms in a unit cell and eight mean fields 
in our theory, (n^^), with r] — A or B indicating the 
sublattice. We use do — 0.334 nm and 27re^(io ~ 3 x 
10-"meV • cm^ 

At Eq — 0, our theory gives a gapped LAF ground 
state at half filled, consistent with the experiment Q, 
and with previous theoretical works by using renormal- 
ization group theoryjs^, H^l and quantum Monte Carlo 
method[34|. First principles calculation indicates that 
such LAF state is stable in the presence of nonlocal 
Coulomb interaction and remote hopping 4l| . In Fig. 
1(a), we illustrate the spin and charge structures of the 
LAF state. The charge distribution is uniform and the 
spins are anti-parallel to each other. There are two- 
fold degeneracy for the spin configurations, related to 
the time reversal. The net spin polarizations in the top 
and bottom layers have an opposite sign, and their sum 
gives a null magnetization. In Fig. 2(a), we show the 
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FIG. 1: (Color online). Schematic illustration of spin (arrows) 
and charge (circles) structures of the LAF state in BLG. A 
and B indicate sublattices, 1 and 2 are layer indices, (a) At 
half filled and E=0, the state is invariant under both com- 
bined transformations T ® I and I ® C, with T, /, and C 
the time reserval, inversion, and particle-hole transformation 
operators, respectively, (b) At half filled and finite E, the 
state is invariant under I ® C. The charge structure breaks 
/, but the spin structure remains invariant under T ® I. (c) 
Charged (electron) graphene in electric field: spin structure 
breaks T ® I and the state is ferrimagnetic. 



LAF energy gap as a function of U . For U « 6.64eV, 
we have a gap of « 2 meV, approximately the gap in 
the BLG measured in the experiment. Note that the esti- 
mated value of U here for the BLG is close to the value of 
U ~ 6.2eV estimated to fit the experimentally observed 
energy gap in the ABC stacking trilayer graphene 0]. 
Therefore, the model seems to gives a comprehensive 
picture about the energy gap of the few layer graphene 
systems 0,111. 

We proceed to discuss the effect of the electric field 
to the half filled BLG. The real space spin and charge 
distributions are illustrated in Fig. 1(b). At half filled, 
an electric field breaks inversion symmetry /, and also 
the particle-hole conjugation symmetry C, but conserves 
I®C. This can be seen from the Hamiltonian for iJp . The 



half filled LAF state in the electric field is invariant under 
I(E)C. Namely, the charge is invariant under I®C , while 
the spin remains invariant under both T ® I and I ® C . 
The system has no net magnetization. In Fig. 2(b), we 
plot the calculated spin dependent energy gap as a func- 
tion of V. As V increases, the spin-down excitation gap 
increases and the spin-up excitation gap becomes nar- 
rowed. At a critical value V — Vc ~ 1.8mV , there is 
a first order transition from the LAF state to the spin 
symmetric normal state. There is a sudden change of the 
gap dX Vc- At y > Vc the gap increases monotonously 
with V . In Fig. 2(c,d), we show the spin polarization 
and the charge transfer at the four distinct lattice sites. 
Our calculations confirm the symmetry analysis we have 
made about the spin and charge structure. There is a 
charge transfer from sublattice A in layer 1 to sublat- 
tice B in layer 2, and a less amount of charge transfer 
from sublattice B in layer 1 to sublattice A in layer 2 
(inter-layer nearest neighboring sites). In Fig. 2(e), we 
show the screening effect and plot E as b, function i?o, 
according to Eq. (4). We remark that Vc ~ \.%mV cor- 
responds to an applied electric field Eq « 15mVnm^^. 
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FIG. 2; (Color online). Bilayer graphene at half filled, (a) 
Energy gap of H in Eqn. (1) as a function of Hubbard U. 
Shown in inset are low energy bands (blue curves). The pa- 
rameters are U — 6.64eV, t = 3.16eV, and t± — 0.381 eV. 
The bands for U = are plotted for comparison (red curves) . 
(b) Spin-resolved energy gap as function of V. (c) spin polar- 
ization on four distinct lattice sites, (d) Charge transfer as 
functions of interlayer electric potential V. In (c) and (d), 
superindices 1(2) and A (B) are for layer indices and sublat- 
tices, respectively, (e) Electric field i? as a function of applied 
electric field -Bo (solid black line), also shown is E v.s. Eo for 
non-interacting case U — (dashed pink line), (f). Energy 
bands at V = 1.4meV. 



This is in good agreement with the external electric field 
at the phase transition point estimated in the transport 
measurement 7]. 

In Fig. 2(f) we shovif the spin-dependent bands near 
the Fermi level for V = lA meV. Both the lower en- 
ergy valence and conduction bands are for spin-up elec- 
trons, and the spin-down electron bands are further away 
from the Fermi level. Because of its semiconducting na- 
ture, all the valence bands are occupied and all the con- 
duction bands are empty at zero temperature, and the 
ground state is spin non-polarized. The spin-polarized 
low energy excitations can be detected in spin-polarized 
transport experiment, and the gap measured by Ve- 
lasco et al. [7| in the presence of external electric field 
should be spin polarized, in agreement with previous 



work 27|. The spin polarization of the excitation gap 
can be detected in transport measurement with ferro- 
magnetic source and drain electrodes. 

The most interesting case is at /i > but small, where 
we have a metallic state in the background of LAF. The 
metallic state in a perpendicular electric field can be a 
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FIG. 3: (Color online), (a) Energy bands of the BLG in an in- 
terlayer electric potential V — 1.4 mV at a charged graphene 
with electron density Sn ~ 2 x lO^^cm"^. Solid lines are for 
spin-up and dashed ones for spin-down. Model parameters 
are the same as in Fig. 2. (b): schematics of probing magne- 
tization by using torque magnetometry experiment. 



half metal. The spin and charge structures of this case 
is illustrated in Fig. 1(c) and there is a net spin- up in 
the BLG. Away from half filled, there is no particle-hole 
symmetry, and the system in the electric field breaks /, 
C, and / ® C, and the spin structure of the LAF state 
breaks T ® I symmetry and gives a half metal with a 
net magnetization. The case for /i < can be obtained 
by particle-hole transformation and will not be discussed 
further. Experimentally, the shift on the chemical po- 
tential can be realized by tuning the gate voltage, which 
may be controlled independently together with the tun- 
ing of the electric field in double gated BLG device. In 
Fig. 3(a), we plot the low energy bands for the BLG 
at the lightly charge density (5n ~ 2 x 10^° cm^^. The 
other parameters are the same as in Fig. 2(f). This 
electron density is equivalent to 2.5 x 10~^ electron per 
carbon site on the BLG. As we can see from the figure, 
the spin-up conduction band is partially filled but the 
spin-down conduction band is completely empty. The 
system is a half metal with a full spin polarization in 
its carriers. The electric field induces time reversal sym- 
metry broken in this case. The surface magnetization 
per area is M = 5n x = 2 x 10~^(7L/iB/2 per nm^, 
with /iG the magnetic moment of graphene atom, and 
/is the Bohr magneton and the Lande g- factor for 
the graphene, which is about 2-2.5. The magnetization 
is tiny, but possibly detectable by using torque magne- 
tometry experiment, as schematically illustrated in Fig. 
3(b). 

We have studied /j, or 6n dependence of the magneti- 
zation and the spin polarization of the BLG away from 
the half filled. As increases, there is a jump in 8n at a 
critical /j,c to separate the LAF and normal states. To il- 
lustrate this, we consider model parameters given in Fig. 
4, which gives a larger LAF gap at half filled, hence a 
more pronounced half metallic state when it is doped. 
The jump in 5n as a function of ^ can be seen clearly in 
Fig. 4(b). The magnetization and spin polarization are 
plotted in Fig. 4(c) and (d), respectively. At 5n = 0, the 
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FIG. 4: (Color online). BLG in an interlayer electric potential 
V = 0.25 eV for model parameters U = 7.2 eV, f = 3.16 eV, 
and t± — 0.5 eV. (a). Energy bands at an electron doping 
(571 ~ 2 X lO^^cm"^. Blue solid lines are for spin-up and 
red dashed ones for spin-down. (b). Electron density 5n vs 
chemical potential. The jump in Sn reflects first order phase 
transition from LAP to normal states, and the system is phase 
separated in the "jumped" density region, (c). Magnetization 
as function of chemical potential (ql = 2; and (d). Site spin 
polarizations as functions of chemical potential. 

spin polarization P_^^ = ~Ps^, and = -^s^- At 
^i > 0,we find Pj^ > |P|^|, and P|^ > |P]^|, which 
results in a net magnetization and is consistent with the 
band picture in Fig. 4(a). 

In conclusion, we have studied the LAF state in the 
BLG by using a mean field theory. At half filled, a weak 
perpendicular electric field leads to spin polarized low en- 
ergy excitations. In slightly charged graphene, the elec- 
tric field can induce a half metallic phase with conduc- 
tion electron or valence hole completely spin polarized, 
and the system is ferrimagnetic. The required strength 
of the electric field and the range of the electron density 
are within the reach of the present experiments in double 
gated BLG. Detecting the nonzero magnetization should 
be a feasible way to identify the half metallic phase in 
the BLG. Such an experiment will also provide direct ev- 
idence for the LAF state. Our results on the half metallic 
phase in the BLG are robust, and are the consequence of 
the symmetry in the LAF state in connection with the 
inversion symmetry broken induced by an electric field 
and particle-hole symmetry broken away from half filled. 
From the band theory point of view, the low energy spin- 
up and spin-down bands of the LAF state in the half filled 
BLG are split by the electric field, which leads to a half 
metal once the carriers are generated in gating. Finally, 
we argue that such half metallic phase should be found 
in the N-layer graphene with rhombohedral stacking or- 
dering, where the LAF state has been proposed to be the 



ground state as we reported recently [i^. The LAF gap 
due to the electron interaction is larger for larger N, and 
saturates at value of about 20 meV around = 9. The 
half metallic phase in multilayer graphene may be more 
robust than in the BLG, but perhaps more challenging 
to realized in experiments. Theoretical work along this 
line is on-going. 
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